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qj Abstract 

<**] , Non-perturbative numerical lattice studies of the Higgs- Yukawa sector of the standard model with exact chiral sym- 
metry are reviewed. In particular, we discuss bounds on the Higgs boson mass at the standard model top quark mass, 
and in the presence of heavy fermions. We present a comprehensive study of the phase structure of the theory at 
weak and very strong values of the Yukawa coupling as well as at non-zero temperature. 
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I. INTRODUCTION 



The Higgs- Yukawa sector of the Standard Model (SM) describes the generation of fermion masses via the non- vanishing 
vacuum expectation value (vev) acquired by the Higgs field which couples through a Yukawa coupling to the fermions. 
The essential element in this picture is that the coupling of the fermions to the Higgs field is chirally invariant which 
leads to the gauge invariant electroweak sector of the SM in the presence of gauge fields. 

There are two couplings in the Higgs- Yukawa sector. They are associate with the Yukawa and the quartic scalar self- 
interaction operators. These couplings are directly related to the fermion and the Higgs boson masses, respectively. 
In the scenario that these masses are large, the corresponding couplings grow strong, and it becomes unclear whether 
the theory can be analysed using perturbation theory or whether non-perturbative methods must be employed. There 
are indeed examples where the applicability of perturbation theory is questionable. The first is the upper Higgs 
boson mass bound which is based on triviality arguments [l| . Here the Higgs boson mass can become large, resulting 
in a strong value of the quartic coupling such that perturbation theory may not work anymore. The second is the 
lower Higgs boson mass bound which is based on vacuum instability arguments [2 5] . Here it is unclear whether this 
instability is not an artefact of perturbation theory applied at large values of the Higgs field such that an expansion 
around the minimum of the effective potential is not justified anymore. 

It is important to stress that both the lower and the upper Higgs boson mass bounds are intrinsically related to the 
cut-off of the theory. Thus, a calculation of the Higgs boson mass bounds can in turn be used to determine the cut-off 
up to which the SM is valid, once the SM Higgs boson mass has been determined. If, for example, the recent result 
for a scalar particle at the Large Hadron Collider (LHC) [1,0] is confirmed as a SM Higgs boson with a mass of about 
125 GeV, the SM could be valid up to very high energies before violating the Higgs boson mass bounds, see Ref. [§[ 
for a recent analysis at next-to-next leading order of perturbation theory. 

Another example where non-perturbative calculations are necessary is the possibility of a heavy fourth fermion gener- 
ation which would lead to a large value of the corresponding Yukawa coupling. Besides these concrete examples, 
it is conceptually very important to study the Higgs- Yukawa sector in a non-perturbative manner since questions such 
as the phase structure of the model, or the spontaneous breaking of the SU(2) ® SU(2) symmetry which underlies 
the Higgs mechanism are of intrinsically non-perturbative nature. 

The need for a non-perturbative investigation of the Higgs- Yukawa sector of the SM has been realised already in the 
early 1990's. A natural choice of a non-perturbative tool is, of course, Euclidean lattice field theory. However, in 
these early studies, the lattice formulations of the Higgs- Yukawa sector were lacking a chirally symmetric form of the 
Yukawa coupling term. The absence of a chirally invariant Yukawa coupling term in the Lagrangian led to severe 
difficulties in studying Higgs- Yukawa model on the lattice, see Refs. [Il| - |l7| and references therein. 

The situation changed however, when it was realised that -based on the Ginsparg- Wilson relation [H|]- there exists a 
consistent formulation of an exact lattice chiral symmetry [l9| , which allows the chiral character of the Higgs- fermion 
coupling structure of the SM to be preserved on the lattice in a conceptually fully controlled manner. This triggered 
a number of lattice investigations of Higgs- Yukawa like models |20l - [28| . 

In this article, we report on the status of the lattice Higgs- Yukawa model using a lattice formulation that obeys an 
exact lattice chiral symmetry as will be explained in Sec.|H] In Sec. IHII we will provide results for the lower and upper 
Higgs boson mass bounds as well as the resonance parameters of the Higgs boson [2^-29]. We also extend the study 
of the Higgs boson mass bounds to the case of a fourth quark generation [3(3]. This calculation will result in rather 
severe constraints on the existence of a fourth fermion generation. 

This article is organised as the following. In Sec. [Hi we describe the setting of our lattice simulations. Section IIIII 
contains results of our work on the Higgs boson mass bounds in the Higgs- Yukawa model. In particular, we have 
investigated the effects of the fermion mass on these bounds. In Sec. lIVl we present our study of the phase structure of 
the model. These include the bulk phase transitions at small values of the bare Yukawa coupling [23j, |24j , as well as in 
the regime of strong- Yukawa coupling [3l[ . We also show results and the status of our work on the finite-temperature 
phase structure in Sec. IIVD1 Finally, we conclude in Sec. [V] 

All statistical errors we quote in this article were obtained with a jackknife or bootstrap analysis, taking possible 
effects of autocorrelations fully into account. Statistical errors of the results presented in Sec. II V CI have also been 
cross-checked using the method in Ref. [HJ ■ 
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II. LATTICE SETTING AND SIMULATION STRATEGY 



The action 



The Euclidean action of the continuum Higgs- Yukawa model containing one doublet of fermions, denoted as and 
b^ c \ and a complex scalar doublet, <p(°', is 



S cont [^ c ),V (c) ,^ (c) ] = / d 4 x|i (d^) 1 (dV c) ) + \ml^ V ^ \ — | y - V 



4 



(c)t,„(c) 



where cp^ — tT2<p^ (r, are the Pauli matrices), 



/(c) 



1-75 



f(c) 

6 (c) 



(c) 



= P+t (c > = [ z ~y :l ) and similar for 6 fi 

In the above equation, mo is the bare mass, Ao labels the bare quartic coupling, and Ht /b denote the bare Yukawa 
couplings. The superscript, (c), in the scalar and spinor fields indicates that these are dimensionful variables defined 
in the continuum. Here we stress that gauge fields are not included in our study, and we perform calculations for only 
one doublet of fermions throughout this work. Moreover, if not stated otherwise, the Yukawa couplings y& and yt, 
are set equal. 

It is straightforward to discretise the pure-scalar component of the above action to obtain 



4 

E 

a=l 



4 att ^E^-E^^+a+E 

X)(JL X 



(2) 



where x is a site on the space-time lattice. The symbol /t denotes the unit vector in the space-time direction /x. The 
mass parameter, mo = amo with a being the lattice spacing, is dimensionless. The real-valued field variables, 
are rendered dimensionless by a proper rescaling with a, and are defined on all lattice sites. These field variables are 
related to the discretised version of the complex scalar doublet, <p^ c \ in Eq. ([!]) through 



(latt) _ 



$2 + jjjjl 
$4 _ ^3 



(3) 



+ A(«-1) 2 



1 - 2A - 8k 



It is convenient to rewrite the scalar action in Eq. ([5} as 

s^ tt = \ ~ 2k E ^x4>x+p, + E [^xK 

with the change of variables, 

* Q = V2^ Q , Ao = \, rhl 
For the fermions we use the action 

S T =E^ [^ + ^ + C^diag(y i ,y 6 )P + +P_diag(y t ,y 6 )0^ a P 



4>x, 



(4) 



(5) 



(6) 



where y t/b = V2K2/ to / &0 , and 



h,2,3 — — * T 1,2,3, $4 — 12x2 j 



(7) 
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where a summation over a is understood. The dimensionless spinor field ip is, 

/ £(latt) \ 

^ = « 3 ( fe (latt) ) > (8) 

with t( latt ) and &( latt ) being the lattice version of and b^ c \ For the fermion kinetic term, we use the overlap 
operator [33-35$ , 

where p is a free, dimensionless parameter, restricted to < p < 2r. The locality properties of the overlap operator 
are optimal for p = 1 in the case of vanishing gauge couplings [Hj], and therefore we set it to this value in this work. 
The operator D w denotes the Wilson Dirac operator defined as 

d w = J2^K-IK% ( 10 ) 

where V£' b,s are the (respectively) forward, backward and symmetrised lattice nearest-neighbour difference operators 
in direction p, and the Wilson parameter r is chosen to be r = 1. The modified chiral projectors are given by: 

p± = l A^, f=-f(i--D<A. (ii) 

This action now obeys an exact global SU(2)x x U(l)y (with Y being the hyper-charge) lattice chiral symmetry with 
the transformations: 

HW + W-^, i> -> i>P+nlulr + rpP-U^, cf> -> u Y (f>n{, ^ ^ n L( j) f ul, (12) 

for any VL L G SU(2) L and U Y € U(l)y. 

B. Implementation 

The actions in Eqs. ((31) and ^ are used in our numerical simulations. We perform calculations on asymmetric 
4-dimensional lattice volumes 

Vi=L 3 s xL u (13) 

where L s and L t are dimensionless spatial and temporal lattice sizes, respectively. In all our zero-temperature 
computations, we choose 

L t = 2L S = 2L, (14) 

with L typically ranging from 8 to 32. We stress that it is essential to perform computations for the Higgs- Yukawa 
models on large volumes. This is because the Goldstone bosons are (almost) massless and induce significant finite-size 
effects proportional to L~ 2 , in contrast to the exponential effects known for a single-particle spectrum and matrix 
elements for theories such as QCD with massive quarks. Figure [T] shows some examples of finite- volume effects that 
are present in quantities investigated in this work. It is clear from these plots that finite volume effects can be very 
large in the calculation of the Higgs boson mass, while they may be mild in other quantities. 
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FIG. 1: Finite volume effects in the magnetisation as defined in Eqs. (|15|l and (|16[) (left), the fermion mass (middle), and the 
Higgs boson mass (right) at a cut-off around 1.5TeV. The data are obtained at infinite bare scalar-quartic coupling, A, and 
fermion masses in the range m/ w 200 — 700GeV. The lattice sizes used are L = L a = 12, 16, 20, 24, 32. We show linear (solid 
lines) and quadratic (dotted lines) fits in 1/L 2 . 



We implement the polynomial Hybrid Monte Carlo (pHMC) algorithm [361 - 4381] . with various improvements (see 
Ref. [3 9| | for a summary) , to perform non-perturbative calculations of the path integral. When compared to simulations 
in QCD using overlap fermions [4(|, it is the absence of gauge fields that makes the application of the overlap operator 
numerically feasible even on large lattices, as it is diagonal in momentum space. 



C. Basic observables 



As described in Sec. Ill Al our simulations are performed using only dimensionless variables in the action. This is 
achieved by rescaling all the dimensionful quantities with appropriate powers of the lattice spacing, a. Therefore, to 
make connection to the real world and to have basic understanding of the spectrum of the theory, it is essential to 
determine the lattice spacing. This is normally carried out by computing the vev of the scalar field, and then setting 
it to the value of 246GeV. Before we describe the details of this procedure, it should be noticed that the scalar vev 
is always zero in a finite system. In principle, one would have to introduce an external source that couples to the 
scalar field and breaks the 0(4) symmetry explicitly, and perform the infinite- volume extrapolation for every quantity 
computed on the lattice, before taking the source to zero. However, this procedure is numerically very demanding, 
and we resort to an alternative method in which we "rotate" the complex scalar doublet in every field configuration, 
such that its ensemble average is given by 



= V / 2rt(m), (15) 



with 



1/2 



x \ a / 



defined on each configuration. It can be shown that the magnetisation, (m), is equivalent to the scalar vev in the 
infinite- volume limit pdT - |43| . 

The rcnormalised scalar vev is given by 

Vr = -£==, (17) 



where Zq is the Goldstone-boson wavefunction renormalisation constant. This renormalisation constant, and the 
Higgs-field wavefunction renormalisation constant Zh, can be extracted from the momentum-space Euclidean prop- 
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agators of the corresponding bosons [H, [28[ , 

t<? Z °l» , (18) 

P + m G/H 

with Oq/h being the Goldstone and Higgs fields, respectively, and all the masses and momenta are in lattice units. 

Through the investigation of the momentum dependence of the Goldstone boson propagator, Z G can be determined. 
This procedure can be improved by performing calculations in one-loop lattice perturbation theory and obtaining the 
propagators to this orderj44(. The lattice spacing, which is related to the inverse of the cut-off scale, A, can now be 
obtained in natural units with, 

. , . 246GeV . , 

o = A _1 , A = . (19) 



The masses of the bosons are given by the pole of the Euclidean propagators in Eq. (|18l) . They can also be extracted 
from the time dependence of the Euclidean correlators with zero spatial momentum |26l l28j , 



:r.y 



AtS>l 



A-g/h exp 



-m G/H L t 



cosh 



m G /H 



At 



(20) 



where A g /h are constants that are proportional to Zq/h. This formula is valid when periodic boundary conditions 
are imposed. Here we stress that this method is applicable only when the ground state is the target single-particle 
state. Therefore, one has to be cautious when studying the Higgs boson, since it may decay into even number of 
Goldstone bosons. The unstable nature of the Higgs boson and the calculation of its resonance parameters will be 
discussed in more detail in Sec. IIII Al 



Finally, to compute the masses of the fermions, we resort to the correlator [26|, |28[ 

C f (At) = ]T ]T (Tr [P^j{t + At, x) ■ ^(t, y)P_ } 



■r.y 



where the trace is over the spinor indices. By studying the time dependence of this correlator, 



Cf(At > 1) oc exp 
the fermion mass can be extracted. 



—mjLt 



cosh 



T~ At 



(21) 



(22) 



III. BOUNDS ON THE HIGGS MASS 



The lattice techniques described in the last section can be applied to the calculation of Higgs boson mass bounds [28, 
[30j . In what follows, we study the model in the broken phase, i.e. where the vev of the scalar field is non-zero. 
The Higgs boson mass is bounded from above by the triviality argument, which reflects the Gaussian nature of the 
fixed point of the theory. This bound is not universal and depends logarithmically on the UV cut-off of the theory. 
Indeed variations in the triviality bound between different lattice regularisations have been observed in the pure 4> 4 
theory [Hj]. 
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There is also an argument from perturbation theory that the Higgs boson mass is bounded from below by a vacuum- 
stability requirement. The picture for the lower bound in perturbation theory arises by examining the effective 
potential. As the fermion fields contribute negatively to the effective potential, they have a destabilising effect. By 
demanding the stability of the theory, this leads then to lower Higgs boson mass bounds. However, it is known that 
the perturbative expansion breaks down for Yukawa couplings near or less than the tree level unitarity bound [46| . 
which is roughly 500 to 600GeV [U H|| . In addition, the perturbative instability occurs at large values of the scalar 
field where an expansion around the minimum of the effective potential may not be trustworthy. Therefore it is 
desirable to have a non-perturbative calculation. 

Although also the lower Higgs boson bound is non-universal, it is expected that it shows a much milder dependence on 
effects of the regularisation employed since a typical ratio A/tuh is of O(I0) for the lower bound, while A/mjj ~ 0.5 
for the upper bound. In the light of the recent discovery of a scalar particle at the LHC, the lower bound becomes 
very interesting: if this scalar particle will turn out to be the Higgs boson, the lower mass bound can be used to 
estimate the breakdown scale of the SM, i.e. the scale where new physics must enter to preserve the stability of the 
theory. 

In this work, we compute the upper and lower bounds of the Higgs boson mass from non-perturbative, direct calcu- 
lations using lattice field theory without rely ing on assumptions such as triviality or vacuum instability. From the 
study of the pure </> 4 theory, it is known [49|-|52j that the Higgs boson mass is a monotonically increasing function 
of the quartic coupling A at fixed lattice spacing. This feature has been demonstrated to be present also in the 
Higgs- Yukawa theory [27| at fixed value of rrif . Therefore, in this work the lower bounds for particular values of m/ 
and A are determined at A = 0, while the upper bounds are obtained at A = oo. 



A. Calculating the Higgs boson mass 



As pointed out in Sec. Ill Cl calculating the mass of the Higgs boson is challenging because of its unstable nature, as it 
decays into even numbers of Goldstone bosons. Extracting the masses and the widths of unstable states in lattice field 
theory is subtle, because the theory is formulated in Euclidean space. It is further complicated by the quantisation 
of spatial momenta in finite volume, since the kinematics may prevent a resonance state from decaying. Therefore, a 
state which is unstable in infinite volume can remain a stable eigenstate in finite volume. 

However, below the inelastic threshold, the infinite-volume phase shift of two-particle scattering can be determined 
via the investigation of finite-size effects in the energy spectrum [53]. Such finite- volume techniques for studying 
scattering states, albeit very challenging to implement in practice, can be used to extract resonance masses and 



widths in Euclidean quantum field theory 54 j 



In this work, we first compute the mass of the Higgs boson by assuming that its width is zero, therefore it is a stable 
particle in finite volume. To check this assumption, we will later use the above-mentioned finite-volume method to 
obtain results of the Higgs boson width, and confirm that the widths is in fact small thus not affecting the results 
assuming a stable Higgs boson. Under the zero width assumption, we extract the Higgs boson mass using the two 
approaches described in Sec. Ill Cl Namely, we study the propagator in Eq. (IT8l) . and the correlator in E g. (l20l . We 
then extract the Higgs boson mass by a fit of the propagator to a perturbation theory inspired formula [26l . |28j and 
by a fit to an exponential form of the correlator of Eq. (|2"0")l . The Higgs boson mass obtained in these two procedures 
are denoted m p H and m c H , respectively. An example of the two methods for determining run is illustrated in Fig. [3] 
We extract the fitted values m p H and m c H which agree within one standard deviation and both fits provide a suitable 
description of the data. The plots in this figure are for m/ = I95GeV. We note that we observe similar agreement 
between m p H and m c H for all our choices of simulation parameters. 

To check the validity of the assumption that the Higgs boson is stable in our work, a calculation of the Higgs boson 
resonance parameters has been performed in Ref. [29j |. Since the finite volume techniques proposed in Refs [53|, [54| are 
only applicable below the inelastic threshold, external sources were introduced which give a mass to the Goldstone 
bosons and break the 0(4) symmetry explicitly. In the calculation the Goldstone boson ener gies were computed at 
non-zero momenta, using the original center of mass frame (53l. IBij as well as a moving frame |55l. [56j. By adjusting 
the values of the external source and the momenta, the Goldstone boson energies were tuned such that 



2E G < m H < 4E G 



(23) 
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The scattering phase shifts from which the resonance parameters were extracted are shown in Fig. [31 along with the 
position of the inelastic thresholds. These phase shifts are used to fit the Breit-Wigner formula to determine the 
resonance mass and width. 
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FIG. 2: Examples of fits to the Higgs momentum space propagator and the Higgs temporal correlation function to obtain 
m v H and m c H , respectively. The results are from a 24 3 x 48 lattice with m/ = 195GeV, A = 1.5TeV. The fitted values are 
m p H — 96.0(4.3)GeV and m c H = 96.4(6.9)GeV where the errors are statistical only and do not reflect the uncertainty in the 
scale determination which, however, affects both values in the same way. 




FIG. 3: Results for the scattering phase shifts at three values of A. From left to right, the plots correspond to > i = 0.001, 10, oo, 
respectively. In each plot, the vertical dotted line indicates the position of the four-Goldstone threshold, above which our 
analysis method is inapplicable. Also, points obtained from both the centre of mass system (c.o.m.) and a system with one 
unit of total momentum (m.f.) are shown. Taken from Ref. [29| . 



The results of the Higgs boson width and mass obtained via the resonance analysis, perturbation theory, using the 
time-slice correlator and employing the momentum space Higgs boson propagator are shown in Tab. HI Here, the top 
quark mass has been set to its physical value. It is clear that the Higgs boson mass determined by the resonance 
study is consistent with that extracted from fits to the momentum space propagator and the temporal correlation 
function. Furthermore, we see that at my = nit, the width of the Higgs resonance is narrow, i.e. at most ~ 10% of 
the resonance mass in all cases. From the results presented in this table, it is demonstrated that it is justifiable to 
assume that the Higgs boson width is zero since it turns out to be very narrow in the resonance analysis such that 
the width has no effect on the mass extraction. 
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A 


A [GEV] 


resonance 


-prcsonancc 
1 H 


-ppcrt 
1 H 


V 

m H 


m c H 


0.01 


883(1) 


0.278(3) 


0.0018(14) 


0.0054(1) 


0.278(2) 


0.274(4) 


1.0 


1503(5) 


0.383(6) 


0.0169(4) 


0.036(8) 


0.386(28) 


0.372(4) 


CO 


1598(2) 


0.403(6) 


0.037(9) 


0.052(2) 


0.405(4) 


0.403(7) 



TABLE I: The results (taken from Ref. 29]) of a study comparing the resonance parameters of the Higgs boson with the results 
of fits to the temporal correlation function and momentum space Higgs boson propagator. Errors are statistical only. Except 
for the cut-off scale, all the results are in lattice units. The fermion mass is set to be the physical top-quark mass. Results 
from three values of the quartic coupling are presented. Also shown are the resonance mass and width from Breit-Wigner fits 
to the scattering cross-section. Finally, a perturbative estimate of the resonance width is included. We note that because of 
some data losses the error on rnP H at A = 1.0 is larger than for the other parameters. 



B. Results of the Higgs boson mass bounds 

We now turn to the results of the Higgs boson mass bound calculations discussed in the previous section. We first 
discuss the results of Ref. [28j], where the upper and lower bounds were computed at several choices of the cut-off 
scale, with the fermion masses at the physical top-quark mass, and also at rrif ~ 676GeV. The main result from 
Ref. [28| is shown in Fig. [4] In the left graph, the situation for a SM top quark mass is shown. The right graph shows 
the situation for a fermion mass of ~ 676GeV. It can be clearly seen that while the upper bound is relatively 
unaffected when using a heavy fermion mass, the lower bound increases substantially. 




500 1000 1500 2000 2500 3000 1500 2000 2500 3000 3500 

Cutoff A [GeV] Cutoff A [GeV] 

FIG. 4: The cut-off dependence of the upper and lower Higgs boson mass bounds for fermion mass at ~ 173GeV (left) and 
~ 676GeV (right). All data have been extrapolated to infinite volume. 



Apart from the cut-off dependence of the bounds at a fixed value of m/, the dependence of the bounds on m/ itself 
has also been examined at a fixed value of the lattice cut-off [57|, the results of which are shown in Fig. [5] (left). We 
clearly observe the increase of the lower bound with increasing m/ in this figure. In particular, Fig. [5] suggests that 
with a Higgs boson mass of ~ 125GeV, the mass of a mass-degenerate fourth generation of quarks is restricted to be 
less than ~ 350GeV. This is clearly already below the bounds from direct experimental searches. 
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FIG. 5: Left: The dependence on the fermion mass of the upper and lower Higgs boson mass bounds, at the cut-off scale 
A = 1.5TeV. Data points from lattice calculations are shown. Results for the lower bound without infinite- volume extrapolation, 
using only 24 3 x 48 lattices, are also shown for comparison. The solid line results from a one-loop calculation of the effective 
potential, as explained in the text. Right: effects of a </> 6 operator with coupling Ae for the lower bound of the Higgs boson 
mass, at various fermion masses and the cut-off scale A = 2TeV. Three values of the coupling constant Ae are plotted. 



In addition to the numerical results, Fig. [5] also contains the estimate of the lower bound from an effective potential 
calculation, which was performed using the same lattice regularisation as in our Monte Carlo simulation. In this 
calculation, the effective potential was computed to one-loop order in the large— Nf limit. Operationally, the one-loop 
calculations were carried out by numerically computing the required momentum-mode summations in a series of finite 
lattice volumes, and then extrapolating to the infinite-volume limit. From this one-loop effective potential, V , the 
Higgs boson mass is determined by solving for the scalar vev, v, and the Higgs boson mass in the gap equations, 

A vm=v=0: -^Vm^ v = m^. (24) 

To compare to the numerically computed lower Higgs boson mass bound, in the effective potential calculation the 
quartic coupling has been set to zero. In addition, the cut-off and the fermion mass were fixed to the same values 
as in the simulations such that a direct comparison is possible. For a standard model top quark mass it has been 
demonstrated in [2(| H3] that the lattice effective potential provides an excellent description for the numerical data 
for the lower Higgs boson mass bound. 

The left panel of Fig. [5] clearly demonstrates that the trend of an increasingly higher value of the lower bound with 
increasing fermion masses, as suggested by the perturbative calculation is realised by the data up to very large values 
of m/ , although the quantitative agreement is better at low m/. Based on this qualitative agreement, we can examine 
the effect of higher-dimensional operators in the effective potential using the same loop and 1 /Nf expansion. To this 
end we include the contribution from the operator Xe<p 6 m the effective potential with Xg the coupling constant. The 
addition of such an operator in the Lagrangian modifies the solution to Eq. (|23|) . and can therefore alter the lower 
bounds on the Higgs boson mass in principle. 

Here we stress that the cut-off cannot be removed in the Higgs- Yukawa model. Furthermore, any perturbative 
expansion in this model is only valid in the regime where the cut-off scale, A = 1/a, is large enough compared to 
low-energy scales such as the Higgs boson mass and the scalar vev. In Ref. [40], it was demonstrated that m/A < 0.5 
(with m being a typical low-energy scale) is enough to ensure the applicability of perturbation theory to the pure 
</> 4 scalar field theory. Here we impose the same condition, but on the value of the scalar field, in our perturbative 
calculation for the effective potential for the Higgs- Yukawa model including the X§4> 6 operator. This results in the 
stability criterion 

J^V(<t>) > Q, 0<O.5, (25) 



where <f> has been properly rescaled to be in lattice units. 
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In the right panel of Fig. [SJ we show the results of our investigation of the lower bounds on the Higgs boson mass, 
using the one-loop effective potential including the contribution from the Xgcf) 6 operator. It is clear that in the regime 
where the perturbative expansion is valid, a wide range of values of A6 lead to qualitatively very similar results. 
Finally, we also point out that exploratory numerical Monte Carlo simulations which include the cf> 6 operator agree 
with the perturbative results for a large range of bare Yukawa couplings [13] ■ 



It is an important task to explore the phase structure of the Higgs- Yukawa model to identify the phase structure of 
the theory and determine the critical coupling constant values where a continuum limit can be performed. In this 
section, we will discuss two aspects concerning the phase structure of the Higgs- Yukawa model considered here. The 
first are the locations of second-order bulk phase transitions in the bare parameter space which can be identified as 
the continuum limits of the lattice theory. For weak values of the bare Yukawa coupling the phase structure has 
been investigated in [23l . |24| and its knowledge was very helpful to identify the simulation parameters for the desired 
physical situation, i.e. a fixed value of the cut-off and the physical values of the fermion masses. Here we remark 
that the bounds on the Higgs boson and fermion masses as presented in Fig. [3] in Sec. IIIIBl are obtained in this weak 
bare Yukawa coupling regime. In this section, we focus now on the large bare Yukawa coupling region and explore 
the phase structure of the theory in this regime of the parameter space. The aim is to investigate, whether the 
phase transitions at large bare Yukawa coupling are governed by the same, Gaussian fixed point as at small Yukawa 
coupling. If we would find deviations from the Gaussian fixed point behaviour, this would open the possibility that 
the renormalised Yukawa coupling can remain strong up to a high cut-off scale which could lead to heavy fermion 
masses and even the existence of bound states. We have therefore been performing simulations at large values of bare 
Yukawa coupling 1 , and the exploratory results will be presented in Sec. IIV CI As a second aspect, we will present 
an investigation of the finite-temperature phase transition in understanding the role of, in particular, heavy fermion 
masses for the electroweak phase transition, especially with respect to questions concerning baryogenesis [60|. 

Before detailing our on-going studies of the bulk and thermal phase transitions of the Higgs- Yukawa model in the 
following two sections, here we describe the general strategy in this work. 

It is natural to use the vev of the scalar field to probe the phase structure. However, a naive computation of this vev 
will always lead to vanishing results in lattice calculations even in the broken phase, because of the finite volume as 
used in the simulations. As discussed in the beginning of Sec. Ill C[ it is appropriate to replace the scalar vev with the 
magnetisation as defined in Eqs. (IT51) and (jTHJ) . 

In order to probe the nature of phase transitions, we have to determine anomalous dimensions of the operators allowed 
by the symmetries. In finite volume, second-order phase transitions are washed out and become cross-overs, and the 
correlation length cannot exceed the size of the system. Therefore, for the study of the phase structure, we resort 
to finite-size scaling techniques. These techniques were developed originally by solving the renormalisation group 
equation (RGE) for finite- volume lattice systems in condensed matter physics [6l| . To draw analogy between field 
theory and statistical mechanics, we also refer to these anomalous dimensions by calling them critical exponents in 
this article, as usually done in statistical mechanics. 

It is challenging to determine the anomalous dimension of the operator corresponding to the Yukawa coupling term, 
because of the presence of fermions and the flavour-changing structure of the operator. We will postpone the discussion 
of this operator for future reports. Here we focus on critical exponents in the scalar sector. To start, we calculate the 
susceptibility, 



IV. STUDY OF THE PHASE STRUCTURE 



A. Purposes and strategy of the study 




(26) 



In Ref. [58, 59], it was demonstrated that in the limit where the bare Yukawa coupling becomes infinity, the Higgs- Yukawa model is 
equivalent to the pure 0(4) scalar model. However, our simulations are performed away from this limit. 
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which is the connected two-point function in the scalar sector. This quantity is proportional to the square of the 
correlation length, £, and diverges at second-order phase transitions in the infinite-volume limit. Solving the RGE for 
this correlator for a finite-size system at fixed cut-off scale (lattice spacing) near a second-order phase transition, one 
obtains the scaling law, 



where g is a universal scaling function, L s is the spatial extent of the lattice, and T c is the critical temperature in 
the infinite- volume limit, which could also be represented by the critical value of a particular coupling. The critical 
exponents, 7 and v are related to the anomalous dimensions of the scalar field and the mass operator, <f> 2 . This scaling 
behaviour is exact near the critical point for space-time dimension, d < 3. Therefore it is an appropriate tool in our 
study of the finite-temperature phase transition. However, in the investigation of the bulk phase structure, we have 
a d = 4 field theory, and the above scaling relation should be modified because of triviality (62-;6(| , if the transition 
is governed by a Gaussian fixed point. These modifications appear as logarithmic corrections in Eq. (|27p . They are 
not included in the analysis presented in this article, but are being considered in our on-going work. 

As will be discussed in the following, the scaling tests and the extraction of anomalous dimensions using Eq. (|27|) arc 
complicated because of the number of free parameters that are involved in the methods for modelling the unknown 
universal function, g. In particular, it is difficult to accurately determine v using this procedure. This complication 
can be reduced by studying Binder's cumulant [67j . 



This quantity is simply the connected four-point function, normalised by the square of the two-point function, in the 
scalar sector. Because of the normalisation, Ql is independent of the critical exponent 7. Furthermore, it is related 
to the renormalised scalar quartic coupling in the infinite- volume limit by a proportionality factor V^/^ 4 [68| . Since 
Binder's cumulant is normalised to be dimensionless, its values computed on different (dimensionless) lattice sizes 
with the same cut-off scale will coincide with each other at the critical point. It is also expected to exhibit milder 
scaling violations resulting from higher-dimensional operators [H, [7(| ■ 

In the next three sections, we discuss details of the investigation of the thermal and bulk phase structures using the 
quantities defined in this section. Errors on all the numerical results in this section are statistical only. 



Before reporting the details of our on-going investigation in the bulk phase structure of the Higgs- Yukawa model in 
the strong- Yukawa regime, we briefly summarise the results obtained in the region of weak- Yukawa coupling [2~ij in 
this section. The order parameters characterising the different phases are the magnetisation defined in Eqs. (1151) and 
(fTB)) , and the staggered magnetisation 




(27) 



B. Bulk phase structure at small Yukawa couplings 




(29) 



The staggered magnetisation is relevant for the breaking of the symmetry, 




(30) 



in the action in Eq. 



In the Higgs- Yukawa model, four phases have been observed: 



1. A symmetric (SYM) phase with (m) = (s) = 0. 
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2. A broken, or ferromagnetic (FM), phase with (m) ^ but (s) = 0. 

3. A staggered-broken, or anti-ferromagnetic (AFM), phase with (m) = but (s) ^ 0. 

4. A ferrimagnetic (PI) phase with (m) ^ and (s) ^ 0. 



Our current knowledge of the phase structure of the Higgs- Yukawa model in the weak- Yukawa regime is summarised in 
Fig. El To make it convenient in comparing results from numerical simulations to a large— Nf analytic calculation |23| . 
we have performed the change of variables, 



UN 



k = k n , A=-i-, $ = WNr$, 



(31) 



in the plots in this figure. The large— Nf calculation was carried out in the Nf — > oo limit while keeping y^, Xn and 
$ fixed. The left panel of Fig. [5] is the result from the large— Nf calculation, and the middle panel is the comparison 
between this calculation and the numerical results from lattice simulations at Nf = 10. The right panel of this figure 
shows the Nf dependence on the critical values of k at the SYM— FM and FM— AFM transitions in our numerical 
calculation, with the Yukawa coupling set to yjv = 0.1. It is observed that the Nt dependence appears to be mild. 
This indicates that the large— Nf analytic calculation may serve as a reasonable, qualitative, guide in choosing the 
simulation parameters for the numerical simulations. Although this analysis has been performed in the weak Yukawa 
coupling region, the good qualitative description makes it possible to also use the large Nf expansion also in the 
strong- Yukawa regime, which was indeed observed in [23| . 
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FIG. 6: The phase structure analysis. The left panel shows analytical predictions for the case of L s = oo, Nf = oo and 
Ajv = 0.1. The black line indicates a first order phase transition, while all other transitions are of second order. The middle 
panel demonstrates a numerical test of the transitions from the SYM to both FM and AFM phases with Nf = 10. The right 
panel displays the Nf dependence in the critical values of k for the SYM— FM and SYM— AFM transitions, at Ajv = 0.1 and 
yN = 1.0. These critical k values are denoted as K^it (> 0) and tt^ rit (< 0), respectively. The squares and circles in the middle 
and right panels of the figures come from direct numerical simulations on the indicated lattice sizes. 



In the weak Yukawa coupling region, we concentrated on the study of the SYM— FM phase transition, which was 
confirmed to be second-order. This allowed us to study physically interesting quantities, such as the Higgs boson 
mass bound presented in Sec. IIII1 near this phase-transition with good control of the cut-off dependence. 



C. Bulk phase transition at large Yukawa couplings 



It is not well understood how the Higgs- Yukawa model at large bare Yukawa couplings differs from that in the weak- 
coupling regime. A first step in a detailed analysis and hence a deeper understanding of the model in this region 
is the investigation of the bulk phase transitions. It can be shown that the Higgs- Yukawa model reduces to a pure 
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scalar non-linear cr-model at infinite bare Yukawa couplings [58|, y59( , and hence becomes trivial at a certain cut-off 
scale. However, it is not clear what happens at large but finite Yukawa couplings. To be able to detect any differences 
from a Gaussian (trivial) theory the critical exponents of the phase transition have to be extracted and compared 
with those of the 0(4) model. If the strong-coupling regime is indeed different from the weak-coupling one and hence 
would be governed by a non-trivial fixed point 2 , it would be very interesting to investigate the possibility of very 
heavy fermions which give rise to a fourth generation, while still maintaining a light Higgs boson in the theory. In 
such a scenario it is unclear, whether an analysis as, e.g. [72[ is applicable and also, whether the Higgs boson mass 
bounds of section III are valid. 

The magnetisation, defined in Eqs. (|15[) and (|16p . can act as an order parameter to identify and determine the order 
of the phase transition. In Fig. [71 the magnetisation for the Higgs- Yukawa model obtained on different lattice volumes 
is shown as a function of y for two n values. In addition, we show the magnetisation as a function of K for the 0(4) 
model. The SYM and FM phases can be clearly distinguished and the phase transition is washed out because of finite 
volume effects as previously discussed. 

The absence of any discontinuities in the magnetisation is strong evidence for a second-order phase transition in all 
three depicted cases. In general, second-order phase transitions are classified through their critical exponents and the 
question arises if these exponents are different in the strong- Yukawa and pure 0(4) models. To answer this question, 
a careful investigation of the susceptibility and Binder's cumulant will be presented in the following. 




FIG. 7: Magnetisation, (m), for the Higgs- Yukawa model at k — 0.06 (left), k — 0.00 (middle) and the pure 0(4) model (right) 
for various volumes. For the 0(4) (m) is plotted as a function of decreasing k to match optically with the Higgs- Yukawa model. 
The absence of discontinuities in (m) is an evidence for a second order phase transition. 



The critical exponents can be calculated by using the finite-size scaling of the susceptibility, Eq. (1261) . The susceptibility 
is shown in Fig. [5] for the Higgs- Yukawa and 0(4) models. This quantity diverges at the critical point in the infinite 
volume limit. Such a divergence in infinite volume is reflected in a bulk finite-size scaling behaviour in lattice 
calculations. As mentioned before in Eq. (f2~T|). the finite-size scaling is predicted by renormalisation group theory, 
with modifications resulting from scaling violation such as that discussed in Ref. (6lj . 



Xm (t,L)-L-^ = g(iLl/^ 



with t = 



t i to 1 



-CL 



-b 



- 1 



(32) 



where C is a phenomenological parameter and b is a shift exponent 61]. This modification comes from the fact that 
the position of the maximum of Xm is volume dependent. From Eq. (1271) the infinite- volume critical temperature can 
be extracted directly. For the 0(4) model we do not observe any shift of the maximum and hence Eq. (|2"T1) is a good 
description of our data in this case. It should be stressed, that the temperature, T, in this section is the control 
parameter. In our work, it is either the Yukawa coupling, y, in the Higgs- Yukawa model or the hopping parameter, 



2 There has been early lattice work on the 3-dimensional Higgs- Yukawa model [7l| , attempting at finding fix points that are different 
from that of the pure scalar field theory. 
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K, in the pure 0(4) model. To extract the critical exponents from the susceptibility, we perform a simultaneous fit of 
all data obtained at all volumes to the partly-empirical formula (73[, 



= A I L- 2 ' u + B t/T( l=co > - C ■ L- b - 1 



-7/2 



(33) 



This formula was also used for a fit to Xm of the 0(4) model, but with the modification of excluding the parameters C 
and b because of the reasons mentioned above. The fit results are summarised in Tab. Inland will be discussed later. 
Notice that there may be logarithmic corrections to the scaling behaviour of the susceptibility because triviality may 
still be present also in the strong- Yukawa model. These corrections should, in principle, be included in Eq. (|33[) 3 . 
This is on-going work, and the result will be presented in a later publication. Therefore, we consider our present 
values of the critical exponents as preliminary and they should be taken with caution. 




FIG. 8: Susceptibility \m at ft = 0.06 (left), n — 0.00 (middle) and the 0(4) model (right) for various volumes. The curves are 
the result of a fit to Eq. (|33|) . The right top boxes in the middle and the right panels show \ m for the largest volumes. For the 
Higgs- Yukawa model a volume-dependent shift of y c towards yi L ~°°' ) can be observed. This shift is not observed in the 0(4) 
model. 







V 


7 


C 


b 


fit interval 


k = 0.06 


18.119(67) 


0.576(28) 


1.038(30) 


4.7(1.6) 


1.95(18) 


17.5, 20.0 


k = 0.00 


16.676(15) 


0.541(22) 


0.996(15) 


10(2) 


2.42(10) 


15.0, 19.0 


0(4) 


0.304268(27) 


0.499(12) 


1.086(19) 


N/A 


N/A 


0.300, 0.308 



TABLE II: Results of a correlated fit to the susceptibility according to Eq. (|33[) where the last column indicates the fit interval. 
The parameter T stands either for y in the Higgs- Yukawa model or for k in the 0(4) model. Since no volume-dependent shift 
can be observed in the 0(4) model for \ m , the parameters C and b have not been fitted here. All quoted errors are statistical 
only. 



It is possible to re-scale the susceptibility according to Eq. (1321) for the Higgs- Yukawa theory, or Eq. (|27[) for the 0(4) 

model, respectively. The fitted parameters extracted from Eq. (l33|) can be used to construct Xm (t,L s ) ■ LJ 1 ^ and 

test its scaling against t- L}J V . This is shown in Fig. [9] Points for all volumes collapse on the same curve in each of the 
three cases shown. This behaviour is typical for second-order phase transitions and hence provides further evidence 
that such a second-order transition happens in the regime of strong Yukawa couplings. 



3 These logarithmic corrections are surely present in the finite-size scaling behaviour of the susceptibility in the pure 0(4) model [62-66]. 
However, our exploratory numerical results show that their inclusion produces minor changes in the results of the critical exponents in 
the 0(4) model. 
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An alternative way of determining critical exponents is via Binder's cumulant, Eq. (|28[). One advantage of this 
quantity over the susceptibility is its milder power-law scaling violation which is given by 

Ql = 9Q L (tL l ' v ) , (34) 

where qq l is a universal function and t is defined in Eq. (|27p . This behaviour can be observed in Fig. [TOl where all 
volumes intersect at the phase transition point in infinite volume where t = 0. Even for the Higgs- Yukawa model no 
shift can be observed and hence the parameters C and b can be completely neglected in the scaling variable. 




FIG. 9: Scaling behaviour of susceptibility at n — 0.06 (left), k — 0.00 (middle) and the 0(4) model (right) for various volumes. 



The value of Binder's cumulant in the broken phase comes from the fact that (m 4 ) w (to 2 ) 2 and hence Ql m 2/3 ^6l\ . 
Our results for Ql at the critical point come close to this value for all setups considered here. Still, Ql obtained in 
the Higgs- Yukawa model differs from the one in the 0(4) model. This may arise from effects of finite renormalisation 
because of the inclusion of fermions. Its implication in the difference of the 0(4) model and the Higgs- Yukawa model 
is under investigation now. Furthermore, it can be demonstrated that for Binder's cumulant, as contrary to the 
susceptibility, there is no logarithmic corrections to the scaling behaviour arising from triviality in the pure 0(4) 
model [rMj . Whether or not such corrections can be present in the Higgs- Yukawa model is being studied now. 




FIG. 10: Binder's Cumulant Ql at K = 0.06 (left), k = 0.00 (middle) and the 0(4) model (right) for various volumes where the 
subscript L indicates the finite volume quantity. Note that the value of Ql at the critical point is different in the Higgs- Yukawa 
and the 0(4) models. 



The basic idea of extracting the critical exponent, v, from Binder's cumulant is the use of the curve collapse of 
Eq. (132|) . If the scaling function gq L is known one will simply minimise [74j 



Rc 



= — \ Ql — 9Ql (tL 1 /^ , (N = total number of data points) 



(35) 
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which would allow to extract v as a direct consequence of the scaling behaviour, 
points, and Rq l is minimal for the correct choice of the parameters v and T C L= °°. 
and systematic errors the function Rq l would become zero. 



The sum is taken over all data 
In the absence of any statistical 



The scaling function <?q l is unknown. However, this can be overcome by the observation that any volume, in the 
following called p, can act as a reference function for the correct choice of parameters, taking thus over the role of 
gQ L . Instead of minimising Eq. (|35|) , we minimise [74| 



Pb = 



■EE E| 

p j^p i,over 



• ft T 1/v \ 
p \~v j J 



1/2 



(36) 



Here, the scaling function is replaced by the interpolating function £ p which is constructed by interpolating the data 

points obtained on volume p to volume j for the values of the scaling variable tijL 1 ^ , with the index i going through 
all data points of volume j. In our case, £ p is computed by picking a point in j and taking the four nearest points 
in p as a basis for a quadratic interpolation. The normalisation factor N oveT is the total number of points used to 
evaluate £ p . The results are summarised in Tab. IIIII and the corresponding curve collapse for Binder's cumulant is 
shown in Fig. [TT] 

In principle, this method could also be used for Xm, but it would be necessary to minimise for five parameters. 
Our investigation shows that this leads to numerical instabilities and the extraction of critical exponents from the 
susceptibility using this method is not possible hitherto. 





rp(L=Oo) 


V 


interval 


k = 0.06 


18.147(24) 


0.550(1) 


17.4, 18.8 


k = 0.00 


16.667(27) 


0.525(6) 


16.0, 17.2 


0(4) 


0.3005(34) 


0.50000(3) 


0.294, 0.314 



TABLE III: Curve collapse results of Binder's cumulant where the last column indicates the interval of the control parameter 
in which the procedure has been used. The parameter T stands either for y in the Higgs- Yukawa model or for k, in the 0(4) 
model. All errors are statistical only. 




FIG. 11: Scaling behaviour of Binder's cumulant at k — 0.06 (left), k = 0.00 (middle) and the 0(4) model (right) for various 
volumes using the parameters listed in table HTT1 



At this point we can claim that we have found a second order phase transition between the SYM and the FM phases 
in the strong Yukawa coupling regime. The absence of discontinuities in (m) and the second-order finite size scaling 
of Xm are strong evidence for such a statement. It is interesting to compare the critical exponents extracted from the 
susceptibility and Binder's cumulant with the ones of the weak- Yukawa model and the 0(4) model. 
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To be able to make a direct comparison of the 0(4) model with the Higgs- Yukawa model, the same strategy has been 
used to compute observables in the pure scalar sector of the theory. In particular, the same analysis techniques have 
been used. The results of the correlated fit to Xm are summarised in Tab. UH The errors quoted there are purely 
statistical. Investigation of the dependence of the results on the fit interval leads to systematic uncertainties which 
are as large as the statistical errors roughly. It is not possible to claim a significant difference in the critical exponents 
between the Higgs- Yukawa model and the 0(4) model from this method so far. 

The curve collapse method, however, can only provide us with information about one critical exponent, namely v. 
The advantage of this method is the significantly smaller statistical error compared to the fit to Xm- However, it 
must be used with care. The scaling behaviour described in Eq. (IMl) is only true close to the critical point. If this 
method is applied at points too far away form the phase transition the result can be affected by scaling-violation 
effects. One possibility to achieve an impression of these effects is the dependence on the interval in which the curve 
collapse method is applied. It was found, that the systematic uncertainty is roughly a factor of five larger than the 
statistical error. However, in the case of k = 0.06 and of the 0(4) model the total error is still a factor of five smaller 
compared to the fitting procedure. In the case of n = 0.00 the total errors are compatible. 

The results of the critical exponent, v, in Tabs. [TT1 and IIIII indicates that the strong- Yukawa model and the 0(4) 
model may belong to different universality classes. However, in the procedure of using Eq. (f3"2j) to determine this 
exponent, the difference of the two models can be as small as two standard deviations. We stress that it is also 
important to investigate the scaling violation as pointed out in Refs. [62l - l66l Ir39l Fzol] . In particular, the observation of 
the multiplicative logarithmic scaling violation is directly related to the triviality of the theory (62t - [66l [75j . Presently, 
we are exploring such analyses and performing computations at additional parameter values. In the near future, we 
will therefore be able to see whether the value of v in the strong bare Yukawa coupling regime is indeed different 
from the one of pure 0(4) model. If we would find a significant difference, then it will be important to investigate 
the strong-coupling regime closer and, in particular, a computation of the spectrum of the Higgs- Yukawa model in 
the strong-coupling region will become most interesting. 



D. Finite-temperature phase transition 



One important subject in the study of the Higgs- Yukawa model is the finite-temperature phase transition. In this 
section we describe the status of our investigation of this transition. We are particularly interested in determining the 
critical temperature where the system undergoes a phase transition from the symmetric phase with vanishing scalar 
vev, v = 0, to the broken phase with non- vanishing v. Further interest lies in the determination of the order of the 
phase transition and the critical exponents. Results reported in this article are obtained at two values of the fermion 
mass, m,f rj 175GeV and m/ « 700GeV. 

Choosing the boundary conditions in the Euclidean temporal direction to be periodic for bosonic and anti-periodic 
for fermion fields, the temperature T on the lattice is given by 

1 A 

T = — = — (37) 
aL t L t 

where L t denotes the dimensionless temporal extent of the lattice. For the study of the finite-temperature phase 
transition, we work at fixed bare Yukawa couplings which lead to the desired fermion masses. Results presented here 
are from lattice simulations performed at A — oo. To vary the temperature, we change the value of k at fixed L t . This 
is equivalent to adjusting the lattice spacing while fixing the number of points in the temporal extent of the lattice 
corresponding then to a change in the temperature. 

Our study shows that the finite-temperature phase transitions in the Higgs- Yukawa model are second-order. The order 
parameter is the magnetisation as defined in Eqs. (|15j) and (|16j) . Since the correlation length is never divergent because 
of finite-volume effects, we resort to finite-size scaling techniques to investigate the second-order finite-temperature 
phase transition in this work. In particular, we analyse the scaling behaviour of the susceptibility of the magnetisation, 
Eq. As in Ref. [73j . we fit the susceptibility according to the partly phenomenologically motivated function 

X m ( K ) = A(i; 2 /" + B +/ _(«- ftc ) 2 ) 7/2 , i/ = 0.68, 7= 1.38, (38) 
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where A, B +/ _, and re c are free fit parameters (B +/ _ are coefficients in the broken and the symmetric phases, 
respectively), and v and 7 are the critical exponents of the three dimensional 0(4) model which are expected to 
characterise the second-order phase transition. Note that we use the fit function of Eq. (|38[) with fixed values of the 
critical exponents only to extract the critical value of re, denoted as re c which in turn leads to the evaluation of the 
critical temperature. This approach is different from that used for the investigation of the strong- Yukawa model as 
described in Sec. lIVCl Since re c depends on the spatial volume, we perform simulations on various spatial lattice sizes 
and perform an infinite volume extrapolation using the formula {D is an unknown constant), 

re c (L) = re c (oo) + D ■ L~ v . (39) 

Having extracted re c in the infinite-volume limit, re c (oo), we can determine the lattice spacing at this re value by 
performing zero-temperature simulations at exactly the same choice of couplings, and using Eq. (|19[) . This then 
allows us to predict the critical temperature, T c , through Eq. (|37l) . In order to estimate the systematic effects in T c 
arising from the uncertainty in re c , we also carry out two additional zero-temperature simulations with re values chosen 
to reflect the error on re c . In this procedure, it is very challenging to maintain a constant Higgs boson mass, since it 
depends significantly on the re value. In the results presented in the following, this is indeed the case. On the other 
hand, we observe that the variation in mu is only ~ 15%, resulting in a similar systematic uncertainty for the critical 
temperature. 



1. Finite-temperature study at physical top quark mass 

As the first step, we investigate the case of a degenerate fermion doublet with the quark mass close to the physical 
top quark mass. To this end we fix the bare Yukawa coupling according to the tree- level estimate of y = m t /v r , which 
has been shown to be a good approximation in this region of couplings )28j . We perform simulations at two different 
temporal extents (L t = 4, 6) for estimating the discretisation effects. In addition, three spatial extents, L s = 16, 20, 24, 
are implemented in order to perform the infinite-volume extrapolation. 

The results of the magnetisation at L t — 4 and 6 are plotted in Fig. [T2T a) and Fig. lT27 dh respectively. It is obvious that 
there is a transition from the symmetric to the broken phase for each choice of L t . The corresponding susceptibilities 
are shown in Figs. [EJb) and IT27 e). The L— dependence of k c {L) is well described by Eq. (f39|) . as can be seen in 
Fig. [T2TD. In Fig. IT2T c). we depict the infinite- volume extrapolation of the cut-off scale, A = 1/a, for the calculations 
using Lt = 4 with three re values around the central value of re c . 

Table IIVI summarises the current status of our finite-temperature study at a fermion mass close to physical top- 
quark mass. Results from two analyses using L t = 4 and 6 are shown. These two choices of the temporal extent 
lead to compatible results of T c , indicating that the discretisation effects may be small. It should be noted that the 
Higgs boson mass turns out to be different in the calculations using these two choices of L t . This can be resolved 
with more tuning and performing new simulations. In particular, we are planning zero-temperature simulations with 
larger spatial extent. This will allow us to have better control of the infinite- volume extrapolation. Nevertheless, our 
preliminary results already show that 

m H /T c ~ 1, (40) 

for mj ~ mf hys . 





re c (oo) 


A(k c ) in GeV 


T c in GeV 


mt in GeV 


mu in GeV 


L t =4 


0.301405(38) 


2070(11) 


517.5(2.8) 


172.1(1.5) 


575(25) 


Lt =6 


0.299561(53) 


3052(64) 


509 (11) 


173.6(2.9) 


500(25) 



TABLE IV: Results for our finite-temperature study performed at the physical top quark mass. 
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FIG. 12: Results of our finite-temperature study at the physical top quark mass with the quartic coupling A = oo. Plots (a) 
and (d): The magnetisation for temporal extents of Lt — 4 and L t = 6. Plots (b) and (e): The corresponding susceptibilities 
with the fit function in Eq. (I38[) . Plot (c): Infinite- volume extrapolation of the cut-off scale at zero-temperature, with the k 
values close to the central value of k c obtained in the Lt — 4 case. Plot (f ) : Infinite- volume extrapolation of k c using Eq. 



Note that for the case of zero temperature L denotes V 14 with V4 = L s Lt and Lt = 2L 



2. Status of finite-temperature study at a quark mass of about 700GeV 



In this section we present the status of our work on the critical temperature in the Higgs- Yukawa model with one 
heavy fermion doublet with a mass of about 700Ge V. We follow the same strategy as in the previous section. Here the 
zero-temperature simulations are still in progress. Thus, the lattice spacings for this calculation are not yet available 
to us. 

Results of the susceptibility, and the infinite- volume extrapolation for k c can be found in Fig. 1131 From Fig. 113r d. 
we notice that the values of k c in the L t = 4 calculation are significantly smaller than that in the L t — 6 analysis. 
From the phase structure presented in Fig. [S] and the value of y ~ 2.8, it is clear that the critical value of k is in the 
FM phase of the zero temperature theory, as expected. Comparing k c for the L t = 4 and L t — 6 calculations, we also 
find that the L t = 4 simulations are carried out closer to the FM— SYM phase boundary, and are thus performed at 
larger values of the cut-off. 
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V. OUTLOOK 



In this article we have provided an overview ol non-perturbative lattice calculations of the Higgs- Yukawa sector of 
the Standard Model and its extension with a fourth fermion generation. The phase diagram of the model has been 
studied and a complex and interesting structure has been revealed. At small values of the bare Yukawa coupling the 
properties of the phase transitions are consistent with the standard model expectation I23l \24 . However, we also 
establish an additional phase transition at very large values of the bare Yukawa coupling [23l 2~il l3l| . This offers 
the very interesting possibility to investigate a strongly interacting Higgs- Yukawa model. We performed a detailed 
study of the properties of the phase transitions at strong bare Yukawa coupling and determined the critical exponents 
characterising the phase transitions through a finite size scaling analysis. Although there are presently indications 
that these critical exponents may differ from the standard model ones, at this stage of our investigations it is too 
early to say that in the strong bare Yukawa coupling region indeed a non-Standard-Model-like phase structure exists. 

As an interesting direction we have also examined the Higgs- Yukawa model at non-zero temperature for fermion 



masses ranging from 175GeV to 700GeV [57|. We find that the transition is always of second order and that the 



critical temperature is higher for increasing fermion mass. 

For a Standard Model top quark mass we have established lower and upper Higgs boson mass bounds as a function of 
the (lattice) cut-off of the theory p6l - l28l | . We also performed a detailed resonance analysis of the Higgs boson which 
confirmed that the Higgs boson mass bounds which assumed a stable Higgs boson are not affected by the resonance 
character of the Higgs boson (2t|. Furthermore, we find that the Higgs boson decay width into massive Goldstone 
bosons is never larger than 10% of the Higgs boson mass and in good agreement with perturbative estimates. As a 
consequence of our lattice study of the lower and upper Higgs boson mass bounds within the Higgs- Yukawa sector at 
a physical value of the top quark mass, we can, in principle, estimate the energy scale at which the standard model 
has to break down. 

We extended the study of the Higgs boson mass bound to a possible fourth generation of quarks considering fermion 
masses up to 700GeV [30(. We found that the upper Higgs boson mass bound shows only a moderate shift by about 
20% at such a fermion mass when compared to the bound for a Standard Model top quark mass. However, the 
lower Higgs boson mass bound is altered significantly and can be as high as 500GcV for a fermion mass of 700GeV. 
We complemented our non-perturbative lattice simulations with a lattice perturbative calculation of the lower Higgs 
boson mass bound from the effective potential. We found a very good agreement with the lattice simulation data. 
This enabled us to test the stability of the lower bound against additions of higher dimensional operators. As a result 
we observed that the lower bound is not affected by including such additional operators. This finding puts severe 
constraints on the fourth generation if the particle with a mass of 125GeV seen at the LHC is the standard model 
Higgs boson. 

Let us discuss the consequences of our lattice study of the Higgs- Yukawa sector of the standard model and its extension 
to a fourth fermion generation, assuming that the particle detected at the LHC @, 0] is a Higgs boson with a mass of 
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125GeV. For the standard model such a Higgs boson mass leads to rather small values of the renormalised quartic and 
Yukawa couplings and it seems therefore that the electroweak sector of the standard model can be described perfectly 
within perturbation theory. Therefore, the perturbative analysis of Ref. [8] provides the result that the energy scale, 
up to which the Standard Model can be valid, is very high. Considering the extension of a fourth fermion generation, 
the lower Higgs boson mass bound together with the phenomenological lower bound of the fourth generation fermion 
mass provides very severe constraints on the existence of the fourth generation. 

As a conclusion, our findings suggest that the electroweak theory of the Standard Model is a perfect description of 
particle interaction up to very high energies as discussed in Ref. [8]. Furthermore, a simple extension of the standard 
model by adding only a fourth fermion generation is most likely not realised. However, as discussed in Ref. [76[ 
the addition of a singlet scalar field could change the situation. As shown in Ref. (76[, the lower Higgs boson mass 
bound can be lowered significantly in the presence of such an additional scalar field. Of course, in Ref. [76| only a 
perturbative calculation has been performed for the scenario of adding such a singlet scalar field and non-perturbative 
calculations, such as the ones presented here, to scrutinise this picture are highly desirable. 

We have demonstrated that with lattice field theory techniques generic strongly interacting Higgs- Yukawa theories can 
be studied in a controlled and accurate way. This became possible through a conceptual breakthrough of formulating 
chiral invariant theories on the lattice together with a much improved understanding of systematic effects such as finite 
size effects or determining resonance parameters. Since in addition the existing computing power of present super 
computers is clearly adequate to perform calculations of Higgs- Yukawa models, lattice computations can contribute 
to our understanding of Higgs- Yukawa models, in particular in the strongly interacting regime. 
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